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ABSTRACT 



We present Giant Metrewave Radio Telescope (GMRT), Hi 21cm observations 
of SBS 0335-052E and SBS 0335-052W, a close pair of dwarf galaxies, which are 
further unusual in being the most metal-poor star-forming galaxies known. We present 
images at several angular resolutions, ranging from ~40 to 4 arcsec. These images 
show that SBS 0335-052 is a strongly interacting system, with a faint diffuse Hi 
bridge seen at low resolution, and elongated tails seen at the higher resolutions. The 
overall morphology suggests that the pair represents a major (as both galaxies have 
similar H I masses) merger of extremely gas-rich galaxies, which is currently past the 
first close encounter. The low-resolution velocity field is dominated by the velocity 
difference between the two galaxies and the velocity gradient along the tidal features. 
However, for SBS 0335-052 W at least, at high angular resolution, one sees a central 
velocity field that could be associated with the spin of the original undisturbed disc. 
The two galaxies have very similar Hi masses, but very different optical properties 
and current star formation rates. A possible reason for this is the differing amounts 
of tidally-induced star formation, because of the different spin orientations of these 
interacting galaxies. The highest angular resolution Hi images show that the ionized 
superbubble, identified by Thuan, Izotov & Lipovetsky (1997), in the Hubble Space 
Telescope (HST) images of SBS 0335-052E, is extended along one of the diffuse tidal 
features, and that there is a high-density H I clump at the other end of the superbubble. 
The star formation in SBS 0335-052E occurs mainly in a group of superstar clusters 
(SSCs) with a clear age gradient; the age decreases as one approaches the dense Hi 
clump. We suggest that this propagating star formation is driven by the superbubble 
expanding into a medium with a tidally-produced density gradient. The high pressures 
associated with the compressed material would also naturally explain why current star 
formation is mainly concentrated in superstar clusters. 
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1 INTRODUCTION 

Low-mass star- forming galaxies with metallicties of Zq/IO 
and lowefl (i.e., 12 + log (O/H) < 7.65; the so called 'eX- 
tremely Metal-Deficient' (XMD) galaxies) are very rare in 
th e local Universe (e.g. , iKunth fc Ostlin 2000t ). According 
to lPustilnik et all (|2003h . they comprise less than 2 per cent 
of all known emission-line galaxies. No more than half a 
dozen (of the tens of thousands of emission-line galaxies) 
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are known to have metallicities near the bottom of the 
metallicity-distribution (Z~Z0/35-Z0/25). In this paper, 
we present GMRT Hi observations of the SBS 0335-052 
system, which contains the lowest known metallicity XMD 
galaxies. SBS 0335-052 consists of a galaxy pair, SBS 0335- 
052E and SBS 0335-052W, with a projected separation 
of ~22 kpc (Pustilnik et al. 2001). Both galaxies, in the 
pair, are extremely metal-poor, and have o xygen abun- 
dance s, 12 + log (O/H) ~7.29 (SBS 0335-052E. fizotov et all 
(|l997l ) and references therein) and 7.12 (SBS 0335-052W, 
Izotov, Thuan & Guseva (2005) and references therein), 
respectively. The galaxies are also peculiar in that they 
have very blue colours in their outer parts (indicative 
of small ages even for the older stellar population, e.g., 
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Table 1. Parameters of the GMRT observations 



Date of observations 


2004 November 28, 29, 




December 10, 11, 12 


Field center R. A. (2000) 


03 h 37 m 44.0 s 


Field center Dec. (2000) 


o / // 

-05 02 40.0 


Central Velocity (km s _1 ) 


4040 


Time on-source (h) 


~20 


Number of channels 


128 


Channel separation (km s" 1 ) 


~3.3 


Flux Calibrators 


3C 48, 3C 147 


Phase Calibrator 


PKS J0323+0534 


Resolution (arcsec 2 ) (root mean 


43 x 39 (1.16) 


square (rms) noise in mjy Bm ) 


21 x 19 (0.85) 




9x9 (0.66) 




6x5 (0.55) 




3x3 (0.51) 



iPapaderos et al l (|l998l ) and iPustilnik. Pramskii fc Kniazevl 
(III)), 

and have the bulk of their baryonic material in the 
form of gas. 

Observations of local, low-metallicity galaxies allow one 
to gain insight into the modes of star formation in primeval 
objects, where the physical conditions (e.g., low metallicity, 
large gas fraction, small gravitational potential, etc.) are 
likely to be similar. Hi observations of local XMD galax- 
ies have been presented by us (e . g., N anqay Radio Tele- 



sco pe data by IPustilnik fc Martini (12007!), and GMRT data 
by Chengalur et al.l 1 2006] ) and lEkta. Chengalur, Pustilnikl 
(2006)). In particular, Ekta, Chengalur, Pustilnik (2008) dis- 
cuss the Hi distribution, kinematics and star formation in 
three of the six most metal-poor star-forming galaxies, viz., 
DDO 68, UGC 772 (SDSS J0113+0052) and SDSS J2104- 
0035. In the case of SBS 0335-052, Very Large Array 
( VLA) Hi 21cm o bserv ations have been presented earlier 
in IPustilnik et alj l|200ll ). Our current GMRT observations 
are both deeper, and also cover a larger range of angular 
resolutions than the earlier VLA ones. Throughout this pa- 
per, we assume a distance of 53.6 Mpc for SBS 0335-052. 
At this distance, 1 arcsec corresponds to a linear distance of 
~260 pc. 



2 OBSERVATIONS AND DATA REDUCTION 

SBS 0335-052 was observed in five different runs, on 2004 
November 28, 29, December 10, 11, 12, at the GMRT, with 
a total bandwidth of 2 MHz, and a channel resolution of 
15.6 kHz (which corresponds to a velocity resolution of 
~3.3 km s _1 ). The total on-source time was ~20 h. The 
flux calibrators, 3C 48 and 3C 147, and the phase calibra- 
tor, PKS J0323+0534 were observed at appropriate inter- 
vals. Each of the data sets was reduced in classic Astronom- 
ical Image Processing System (AIPS). The data reduction 
steps included flagging of the bad visibility points, calibrat- 
ing for phase and bandpass shape. The visibility data from 
the different runs were combined into one data set using the 
task 'DBCON'. Since the differential Doppler shift between 
the different runs is small compared to the channel width, 
no Doppler correction was applied. After continuum sub- 
traction, each channel was imaged using task 'IMAGR'. The 
parameters of the GMRT observations are given in Table [1] 



3 H i MORPHOLOGY AND KINEMATICS 

The integrated Hi intensity map of SBS 0335-052 system, 
at an angular resolution of ~40 arcsec, is shown in Fig.QJA). 
At this resolution, the two galaxies appear to be embedded 
in a common Hi envelope. The earlier VLA observations 
(Pustilnik et al. 2001) also show the galaxies to be in a 
common envelope, but the total H I extent seen in our map 
(225 x 105 arcsec 2 , or 58.5 x 27.3 kpc 2 , at an Hi column 
density level of 1.0 x 10 19 atoms cm -2 ) is smaller than that 
deduced from the VLA observations (65 x 22 kpc 2 , after 
rescaling to the distance adopted here). Our Hi maps are, 
at least, a factor of 1.5 times more sensitive than those from 
earlier VLA observations. The rms noise quoted by the lat- 
ter is 1.0 mjy beam -1 per channel, for a channel width of 
5.3 km s _1 , at an angular resolution of ~20 x 15 arcsec 2 . 
Our maps, at slightly coarser resolution (~21 x 19 arcsec 2 ), 
have an rms noise of 0.85 mjy beam - per channel, for a 
channel width of 3.3 km s _1 . The outermost contour, in 
Fig IH A), is a factor of ~5 times lower than the sensitivity 
level in the VLA data. The difference in linear sizes seen at 
the GMRT and VLA is largely because a tentative feature 
(to the north-east of the main emission) seen in the VLA 
image has not been detected at the GMRT. 

The GMRT Hi velocity field, at ~40 arcsec, is shown in 
Fig.[TjB). There is an overall velocity gradient in northeast- 
southwest (NE-SW) direction. The contours are crowded in 
the region between the two galaxies. A contour map (at the 
same angular resolution) of every alternate channel, in the 
data cube, is shown in Fig. [2] H I emission, peaked at the 
location of each of the two galaxies, and an Hi 'bridge', 
connecting the two galaxies (at intermediate velocities) can 
be seen. 

Fig. [T] also shows the Hi emission at two higher angular 
resolutions, viz., ~ 20 arcsec (panels C, D) and ~9 arc- 
sec (panels E, F, G). In the ~20 arcsec image, the 'bridge' 
emission appears to be completely resolved out. To check if 
the 'bridge' emission seen in panel A is entirely due to beam 
smearing, we smoothed the ~20 arcsec map to 40 arcsec res- 
olution, and compared it with the map shown in panel A. 
The smoothed map does not recover all the flux seen in the 
'bridge' region in panel A, indicating that the emission seen 
is not entirely due to beam smearing. At ~20 arcsec reso- 
lution, one can clearly distinguish the Hi associated with 
each galaxy. Both galaxies show extended elongated fea- 
tures, reminiscent of 'tidal tails'. The velocity field of the 
SBS 0335-052 system, at a resolution of ~21 x 19 arcsec 2 
in Fig. HID), shows complex patterns. Both galaxies have 
disturbed fields. In the case of SBS 0335-052E, there is a 
clear NE-SW velocity gradient, but the velocity contours 
show substantial deviations from those expected from a ro- 
tating disc. In the case of SBS 0335-052W, there is an overall 
NE-SW velocity gradient along the 'tidal tails', but in the 
central part of the galaxy, the velocity gradient is approxi- 
mately from south to north. 

At a still higher angular resolution (viz., ~9 arcsec, 
panel E), one can see more clearly that the gas distribution 
of each galaxy consists of a central peak coincident with the 
optical emission, along with elongated 'tidal tail' like fea- 
tures. In the case of SBS 0335-052W, one can also see that 
the central Hi body is elongated in the north-north-west- 
south-south-east (NNW-SSE) direction (see also Fig. [5]). 
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Figure 1. (A.) Integrated Hi emission map of SBS 0335—052 system, in contours, at a resolution of ~43 X 39 arcsec 2 , overlaid on the 
B-band Digitized Sky Survey-II (DSS-II) image (grey-scale, arbitrary units). The contour levels are at Hi column densities of ~0.14, 
0.31, 0.69, 1.52, 3.34 X 10 20 atoms cm -2 . (B.) The Hi intensity-weighted velocity field (contours and grey-scales) at a resolution of 
~43 X 39 arcsec 2 . The velocity contours are from 3997 to 4074 km s~ 1 , in steps of of 3.3 km s — 1 . (C.) Same as (A), excepting that 
the Hi image is at an angular resolution of ~21 X 19 arcsec 2 . The contour levels are at Hi column densities of 0.30, 0.67, 1.50, 3.33 
and 7.40 X 10 20 atoms cm . (D.) Same as (B), excepting that the Hi velocity field is at an angular resolution of ~21 X 19 arcsec 
resolution. The velocity contours range from 3997 to 4080 km s _1 , in steps of 3.3 km s _1 . (E.) Same as (A), excepting that the Hi 
image is at an angular resolution of ~9 arcsec. The contour levels are at Hi column densities of 0.75, 1.39, 2.57, 4.75, 8.78, 16.24 and 
18.81 X 10 20 atoms cm . (F,G.) Same as (B), excepting that the Hi velocity field is at an angular resolution of ~9 arcsec. In F, 
the velocity contours range from 4037 to 4060 km s" 1 (SBS 0335-052E), and in G, from 4013 to 4036 km s" 1 (SBS 0334-052W). The 
contour spacing is 3.3 km s _1 . 
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Figure 3. The spectra of SBS 0335-052E and SBS 0335- 
052W (left) are overlaid upon the spectrum of the whole sys- 
tem. These are obtained from data at an angular resolution of 
~43 X 39 arcsec 2 . 



Figure 2. The channel maps of SBS 0335-052 system, in 
contours, at a resolution of ^43 X 39 arcsec 2 . The contours 
are at -3, 3, 4.2, 6.0 times the rms noise in a single channel 
(1.16 mjy beam -1 ). Every alternate channel is shown. The B- 
band DSS-II image is shown in grey-scale, in arbitrary units. The 
velocity, for each shown channel, is given in corresponding panel. 



The Hi velocity field, at this resolution (panel G), shows 
that the velocity gradient in SBS 0335-052W is also ap- 
proximately aligned with the major axis of the central H I 
distribution. This suggests that at this resolution one is be- 
ginning to see the gas associated with the disc of this galaxy. 
In the case of SBS 0335-052E, the dominant velocity gradi- 
ent still seems to be along the tidal tail. 

To summarize, our multi-resolution images show 
SBS 0335-052 to be strongly interacting, with a diffuse H I 
'bridge' seen at low resolution, and elongated 'tidal tails' 
seen at the higher resolutions. The large-scale velocity field 
is largely dominated by the velocity of the tidal features, 
however, for SBS0335-052W at least, at the highest angular 
resolution, one sees a central velocity gradient that one could 
associate with the spin of the original undisturbed disc. 

The spectra of SBS 0335-052E, SBS 0335-052W 
and the whole system (obtained from ~40 arcsec data) 
are shown in Fig. [3] The total integrated Hi flux is 
1.46 ± 0.15 Jy km s" 1 , whi ch matches, within error 
bars, with the value given in iThuan et all (|l999l ), viz. , 
1.28 ± 0.22 Jy km s" 1 . The VLA observations (Pustilnik et 
al. 2001) obtained a higher value of 2.46 ± 0.18 Jy km s -1 . 
The reason for this discrepancy (which is at the ~4cr level) 
is unclear. We have confirmed the amplitude calibration of 
the GMRT maps, by comparing the fluxes of background 
continuum sources in our map with their NRAO VLA Sky 
Survey (NVSS) fluxes. The Hi fluxes measured for eastern 
and western components are 0.61 and 0.86 Jy km s" 1 (42 
and 58 per cent of the total flux), respectively, correspond- 
ing to Hi masses of 4.2 and 5.8 x 10 8 Mq, and an Hi mass- 
ratio of 1:1.4. The central velocities of SBS 0335-052E and 
SBS 0335-052W are 4053.6 ± 1.7, 4014.7 ± 1.7 km s _1 , re- 
spectively, and the velocity widths, at 50 per cent levels, are 



Table 2. Main parameters of the observed galaxies. 



Parameter 


SBS 0335-052E 


SBS 0335-052W 


V hc i (km s" 1 ) 


4014.7 


4053.6 


12+log(0/H) 


7.29 


7.12 


Hi flux (Jy km s" 1 ) 


0.61 


0.86 


M H i (10 8 M ) 


4.2 


5.8 


M s tar/(M st ar+M gas ) 


0.07 


0.016 


M dyn (10 9 Mq) 


7.0 


7.9 



50.8 ± 3.3 and 47.4 ± 3.3 km s _1 . For the entire system, the 
central velocity is 4031.5 ± 1.7 km s , the velocity range 
over which emission is detected is 73 km s _1 , while the total 
linear extent is ~58.5 kpc. Assuming the system is bound, 
its indicative dynamical mass can be estimated as: 

M ind = 2.3 x 10 5 x R kpc x V 2 km M Q = -9.0 x 10 9 M . 

The same formula gives indicative dynamical masses of 
6.0 x 10 9 M and 7.2 x 1O 9 M for the east (E) and west 
(W) galaxy, respectively (where we have assumed sizes of 
10 and 13.6 kpc, and velocity widths of 51 and 48 km s~\ 
respectively). The sum of the indicative dynamical masses 
of the two galaxies exceeds the total indicative dynamical 
mass, because the various inclination angles involved have 
not been correctly accounted for. 

The indicative dynamical mass is ~10 times larger than 
the total gas mass. Using the stellar masses from Pustilnik et 
al. (2004), the M star /(M s tar + M gas ) ratio is 0.07 and 0.016 
for the E and W galaxy, respectively. The various measured 
parameters for the two galaxies are summarized in Tabled 

4 STAR FORMATION IN SBS 0335-052 

Current star formation in SBS 0335-052E is largely 

confined to several superstar clusters ( SSCs), locate d 

in the central region of the galaxy l|Thuan et al. 19971 ). 
The ages of the SSCs show a systematic variation 
with position, with the oldest SSCs being to the 
north a nd the youngest clusters being at the south- 
ern tip l|Thuan et al. 19971 ; iReines. Johnson fc Hunt 20081 ; 
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iThompson et al. 20091 ). Reines et al. (2008) compute an age 
of ~15 Myr for the oldest SSCs and <3 Myr for the youngest 
ones, and also determine that the star formation is propa- 
gating at a rate of ~35 km s _1 . Thuan et al. (1997) also draw 
attention to a large (~380 pc) superbubble, which they iden- 
tify as being a supernova blown cavity. Fig. 2Jleft) shows the 
Hi emission (contours) at resolution of ~7 arcsec (1.8 kpc) 
overlaid on an HST Advanced Camera for Surveys (ACS) 
image. As can be seen, the Hi contours in the north-west 
are aligned with the direction of the superbubble. The asym- 
metry of ionized gas, r elative to the central sta rburst, is also 
seen on th e data from | Papaderos et al.l l|l998h . Pustilnik et 
al. (2004). Ilzotov et alj (120061 ). and recent Fabry- Perot in- 
terferometry to study Ha-kinematics (Pustilnik et al., in 
preparation). Note also that the north-west extensions of 
the H I contours are essentially due to material being pulled 
out to form the 'bridge' between SBS 0335-052E and W 
(see the lower resolution images, Fig. [T|. At a still higher 
resolution of ~4 arcsec (1 kpc, Fig. Upright)), one can see 
that the detected Hi emission is confined to a region to 
the south-east of the southern most SSC. The sensitivity of 
the Hi observations decreases with increasing spatial reso- 
lution, so the high resolution images detect only the highest 
column density gas. As can be seen from Fig.|3Jleft), the su- 
perbubble region corresponds to gas with column densities 
<9 x 10 20 atoms cm , while the density in the south-east 
ridge (Fig. fright)) reaches up to ~21 x 10 20 atoms cm -2 . 
The superbubble seems to have expanded preferentially into 
the lower density gas. The latter is produced when the 
'bridge' material is being pulled out. The propagation of 
star formation towards the south-east could be due to the 
propagation of the superbubble into the relatively denser 
gas there. 

In this context, it is interesting to note that the expan- 
sion velocity that Thuan et al. (1997) derive for the super- 
bubble is ~17 km s~ , which is within a factor of ~2 of 
that derived for the propagation rate of the star formation. 
Given the large uncertainties in both these rates, agreement 
to within a factor of 2 is quite reasonable. The north-to- 
south propagation of the star formation, hence, appears to 
be attributable to the expansion of a superbubble into a 
medium in which there is an overall north-south density gra- 
dient because of the tidal interaction. Such a scenario would 
also provide a natural explanation for the dominant cluster 
mode of star formation observed in SBS 0335-052E. High 
pressures are expected in gas that is being compress ed by 
an expanding supernova shell, and lElmegreerj |2004l ) sug- 
gests that high pressure is conducive to the cluster mode of 
star formation. A similar pattern of sequential star forma- 
tion and a dominant cluster mode for the star formation has 
recently been found associ ated with an expandin g supernova 
shell in the outer galaxy (|Kobavashi et al1l2008l ). 

Given that SBS 0335-052E and W are among the low- 
est metallicity galaxies known, it is interesting to see if 
the star formation occurs at a similar gas thres hold den- 
sity a s for more metal-rich galaxies. For example, ISkillmanl 
(|l987i ) suggested that star formation in dwarf galaxies oc- 
curs only when the gas column-density crosses a threshold 
value of ~10 21 atoms cm -2 ; this in turn could be related 
to a threshold amount of dust shielding required for the 
production of molecular gas. In this case, one would expect 
that the threshold column-density increases with decreasing 









1 1 1 1 1 1 


8 40 






1 




1 


U 1 










:0 


i i -i 




1 1 





Figure 4. (Left.) The integrated Hi intensity map of SBS 0335- 
052E, at an angular resolution of ~7 X 6 arcsec 2 , in contours, 
overlaid on an HST ACS continuum image (grey-scale). The 
contours are at Hi column densities of 1.3, 2.3, 4.0, 6.9, 11.9, 
20.6 X 10 20 atoms cm . The grey-scale is in arbitrary units. 
(Right.) Same as on the left, except that the resolution of the 
Hi image is ~4 X 3 arcsec 2 , and the contours are at Hi column 
densities of 8.7, 13.2, 20.1 X 10 20 atoms cm" 2 . 





Figure 5. (Left.) The integrated Hi intensity map of SBS 0335- 
052W, at an angular resolution of ~7 X 6 arcsec 2 , in con- 
tours, overlaid on an DSS-II B-band image (grey-scale). The con- 
tours are at Hi column densities of 3.6, 5.3, 7.7, 11.4, 16.7 and 
24.6 X 10 20 atoms cm . The grey-scale is in arbitrary units. 
(Right.) Same as on the left, except that the resolution of the 
Hi image is ^2.9 arcsec, and the contours are at Hi column den- 
sities of 1.6, 2.3, 3.3 and 4.7 X 10 21 atoms cm -2 . 



gas-phase metallicity. Schavej l|2004f ) (see also ISchave 20081 ; 
ISchave fc Dalla Vecchia 2 008 ) suggests a different model, 
one in which star formation is related to the formation of 
a cold phase of the interstellar medium. From his fitting 
formulae for the threshold column density, the threshold 
density increases with increasing gas fraction and decreas- 
ing metallicity (all other factors being equal). In this model 
too, gas-rich XMD galaxies should have a higher threshold 
for star for mation than spirals . This issue was also explored, 
in detail, in lEkta etT al. (2008), who found no conclusive ev- 
idence for an increase in threshold density with metallicity, 
albeit in a relatively small sample of XMD galaxies. 

The H I column densities at the positions of the bright- 
est star-forming regions in SBS 0335-052E and SBS 0335- 
052W are ^2.0 and 5.4 x 10 21 atoms cm -2 , respectively. 
The measurements were made at a resolution of ~3.4 and 
2.75 arcsec (880 and 715 pc), respectively. The signal-to- 



6 Ekta et al. 



noise ratio for SBS 0335-052E was too low for the col- 
umn density to be measured at the higher resolution used 
for SBS 0335-052W. In SBS 0335-052E, the Hi peak is 
slightly (i.e., ~1.5 arcsec, or 0.4 kpc), offset from the bright- 
est star-forming region (as was already noticed from the 
VLA data (Pustilnik et al. 2001)). Nevertheless, it lies 
within a half beam-width of the latter. The western H I 
peak is also slightly (~1 arcsec, or 0.26 kpc) offset from 
the brightest star-forming region in SBS 0335-052W. Note 
that the column densities have not been corrected for in- 
clination, because, as discussed above, it is difficult to es- 
timate the inclination angles of the Hi discs. The Hi col- 
umn density of SBS 0335-052 E, as calculated from its 
HST spectral observations by iThuan fc Izotovl (|l997f ). is 
7.0 ± 0.5 x 10 21 atoms cm -2 . The higher value that they ob- 
tain is presumably a result of the much higher angular reso- 
lution of the HST data. These peak H I column densities, are 
amongst the largest values found in ga s -rich dwarf or XMD 
galaxies (see, e.g. , iTavlor et al. 1994 ; Ivan Zee et al. 19981 ; 
iBegum et al. 20061 ; Ekta et al. 2006; Ekta et al. 2008). This 
is despite the fact that all the above mentioned studies have 
better linear resolution than that of our observations, except 
for a few galaxies in Taylor et al. (1994). Given the uncer- 
tain inclination correction however, it is difficult to gauge 
the significance of this higher observed H I column density. 



5 DISCUSSION AND SUMMARY 

SBS 0335-052 is the lowest metallicity star-forming system 
known. The system is further peculiar in being an interacting 
pair of dwarf galaxies. In the case of SBS 0335-052E, the 
metallicity is low not just in the ionized Hn regions, but 
also in the surrounding Hi gas. From FUSE observations, 
IThuan. Lecavelier des Etangs fc Izotovl (|2005l ) find that the 
metallicity in the H I envelope is similar to that in the central 
star-forming regions. 

Our H I data shows that the pair is in an advanced stage 
of interaction, with the gas being pulled out to form a bridge, 
as well as into elongated tidal tails. As discussed in Sec- 
tion [3] if the dynamic-to-baryonic mass ratio is >10, the 
system is bound, and will eventually merge. The relatively 
short tails, and the faintness of the 'bridge' then suggest that 
the pair is probably past its first close encounter, but has 
yet to fall towards each o ther for the second tim e (see, e.g., 
iToomre fc Toomrdll972T ). Idi Matteo et all IMPOST ) show that 
starbursts can be triggered in extremely gas-rich galaxies, 
which are undergoing their first close encounter, depending 
on the details of the galaxy structure and the orbit geometry. 
The starburst onset could range from being just before the 
first close encounter (e.g., for very close, planar encounters) 
to 0.1-0.4 Gyr after the encounter. The relatively broad tidal 
features (as opposed to the narrow long features produced 
by direct encounters) also suggest that neithe r of the galax- 
ies is undergoing a direct encounter (see, e.g.. lHoward et ah! 
1993). Indeed from Fig.[U one can see that if the elongated 
structure seen in the central regions of SBS 0335-052W re- 
flects the underlying gas disc, the galaxy is undergoing a 
highly inclined encounter. 

Interestingly, all of the six lowest-metallicity star- 
forming XMD galaxies known are either currently under- 
going tidal interactions, or appear to be the remnants of a 



merger of extremely gas-rich progenitors l|Ekta et al.ll2008n . 
While interactions have long been regarded as triggers for 
starbursts, numerical simulations show that, in fact, only a 
small fraction of all binary interactions lead to a reasonably 
strong starburst (jdi Matteo et al J 12008). The enhancement 
in the star formation rate during merger depends on a num- 
ber of factors, but in general (in recent models) the inter- 
nal structure (e.g., presence of a bulge, gas-mass fraction, 
disc-stability factor, etc.) appear to play a more crucial role 
in determining the strength and duration of the starburst 
than the interaction geometry (|Mihos fc Hernquistl 1 19961 ; 
di Matte o et al. 1200811 . In particular, for gas-rich galaxies, 
di Matteo et al.l (|2008h find that the amplitude of the star- 
burst, which occurs just after the first close passage, is gen- 
erally stronger than for gas-poor galaxies. 

In the current situation, the galaxies concerned are both 
extremely gas-rich dwarfs, and are likely to have similar in- 
ternal structure. While numerical simulations of such sys- 
tems are very limit e d, for somewhat less gas-rich systems, 
iMihos fc Hernquistl (|l996l ) show that coplanar orbits lead 
to somewhat stronger starbursts than inclined orbits. As 
discussed above, SBS 0335-052W appears to be undergo- 
ing an inclined encounter, which could perhaps be the rea- 
son for its relatively lower star formation rate. As far as 
simulations of interac t ing, g as-rich galaxies are concerned, 
ISpringel fc Hernquistl |k)05) show that a direct encounter 
between two giant, extremely gas-rich discs can lead to the 
formation of spiral-galaxy-like gas-rich remnant. The star 
formation rate is enhanced both during the first close pas- 
sage as well as during the final merger. In general, gas-rich 
mergers tend to show larger enhancements in the instan- 
taneous star formation rate than gas-poor ones, both be- 
cause there is more raw material available for star forma- 
tion, and also because of the inherent instability of thin gas 
discs (Idi Matteo et al-lbOOST ). 

We used th e rec ently released Identikit tool 
( Bar nes fc Hibbardl [2009) to model the interaction in 
the SBS 0335-052 system. Since Identikit is based on 
collisionless N-body simulations, we compare only the 
morphology of the outer regions and the velocity-field 
projections. The best possible match, that we could find, 
is shown in Fig. [6] along with the data projections in the 
same co-ordinate system. While the match is not perfect 
(given the limitations of the model, a perfect match is 
perhaps not to be expected) the parameters from the match 
could be taken as indicative. In the displayed model, the 
galaxies undergo a close (pericenter distance of ~0.55 kpc) 
encounter, which is somewhat retrograde encounter for 
SBS0335-052E, and close to polar for SBS 0335-052W. 
The sky plane is inclined ~45 degrees with respect to 
the orbit plane. The observed snapshot corresponds to 
shortly (~0.1 Gyr) after the first perigalacticon. From 
optical colours, the major enhanced star formation in both 
components is estimated to have started 0.1-0.4 Gyr ago. 
It is interesting to note that this is consistent with the 
star formation acitivity starting just before the pericenter 
passage, as predicte d for some cases of ver y close encounters 
in gas-rich mergers (|di Matteo et al.ll2008l ). A more detailed 
simulation of this unique pair would be highly desirable. 

While the interaction in SBS 0335-052 is likely to have 
triggered the current starburst, it almost certainly plays 
an important role in the long-term evolution of the sys- 
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Figure 6. Left: Projections of the Hi data cube, made for comparison with the Identikit model. The panels (clockwise and starting 
from the upper left) are (1) The sky (X, Y) plane projection of the Hi in SBS 0335—052 system, at an angular resolution of ~40 arcsec. 
Note that the RA ('X') axis is flipped with respect to Fig. [T] so the E galaxy is to the right. (2) The V, Y plane projection (where 'V is 
the velocity). (3) The X, V plane projection. Right: Identikit web interface results for the closest fit model that we could make for our 
data. The panels are in the same order as for the observed Hi data. The lower right panel is the X, Z plane projection, where Z is the 
depth along the line of sight (for which there is no counterpart in the left panel, since the Z co-ordinate is not an observable). Blue dots 
correspond to the first galaxy, and pink ones correspond to the second galaxy. The input parameters used in the model are pericenter 
distance p = 0.125 (1.5 h), inclinations of galaxy spins relative to the orbital plane ii=130, 12=70 degrees, azimuthal orientations cji=340, 
1^2=90 degrees, the angles between the sky plane and the orbital plane, 0x=45 and <9y=45 degrees. The rotation angle around Z-axis 
is 310 degrees. The snapshot is for a time t ~0.5 (or ~70 Myr, in physical units), after the first close passage. The conversion from 
dimensionless units to physical units was done assuming that the galaxies in SBS 0335-052 system have scale-length, h ~0.37 kpc and, 
V rot ~40 km s -1 . 



tern. iBekkil I 20081 ) presents results from numerical simu- 
lations of gas-rich dwarf galaxies, and shows that in the 
case where there is no stellar bulge, and the gas is ini- 
tially in a cold disc, the merger remnant resembles a mod- 
erately gas-rich BCD (blue compact dwarf) galaxy. The 
central compact emission comes from a starburst triggered 
by the inflow of gas from the outskirts of the initial cold 
disc. High gas pressures in the star- forming regions also 
lead to the formation of massive star clusters (as seen for 
SBS 0335-052E). As mentioned above, the role of high pres- 
sures in forming massive star clusters was al s o disc ussed 
earlier by Elmcgrccn, Kaufman & Thomasson ( 1993|) and 
lElmegreenl i|2004h . Further. ldi Matteo et all (| 20081 ) find that 
the strength of the central starburst found in merging gas- 
rich galaxies depends critically on the assumed disc stabil- 
ity. Simulations in which gas dissipation has been decreased 
to prevent excessive disc fragmentation for isolated galaxies 
(e.g., their particle-mesh - sticky-particles (PM-SM) simula- 
tion) show very little enhancement of star formation during 
the merger. Further, even in their simulations which do show 
an enhancement in the star formation rate, the major cause 
of the enhancement is induced disc instabilities, and not in- 
flow of gas. While a significant radial metallicity gradient is 
known to exist in large spiral galaxies, such strong gradi- 
ents have not been found in dwarf galaxies, and hence, the 
role that infalling gas plays in lowering the metallicity of the 
central regions may not be important. 

Inflow of metal-poor gas is one possible cause for siz- 
able metal deficiency in central star-forming regions of some 
dwarf galaxies. Other possible causes are (a) preferential 
loss of metal-enriched gas, i.e., escape of supernova metal- 
enriched material from previous starbursts, (b) infall of pris- 



tine 'gas clouds' from the intergalactic medium, (c) very slow 
evolution of locally very stable gas discs (extreme low sur- 
face brightness (LSB) galaxies), which fortuitously (e.g., due 
to residing in a void) escaped strong external disturbances 
and related additional star formation, and (d) a genuinely 
young galaxy which has not undergone much star formation 
(i.e., due to a too small an age). 

Numerical simulations (e.g., Fra gile et alj [2004) show 
that sufficiently strong starbursts, in sufficiently small galax- 
ies, could indeed lead to significant outflow of metal-enriched 
gas. A high spatial and velocity resolution study of the ion- 
ized gas around the starburst region in SBS 0335-052E 
has i ndeed identified a high velocity (~50 km s _1 ) out- 
flow |lzotov et al l l2006h . The final fate of this outflowing 
gas (i.e., whether it escapes from the galaxy or falls back 
to form a 'galactic fountain') depends on several parame- 
ters, including the mass of the parent galaxy, the presence 
of a hot-gas halo from the previous star formation activity, 
etc.. In addition, the positive feedback from massive, com- 
pact superstar clusters could also inhibit outflow of g as (e.g., 
iTenorio-Tagle et alj|200l 120071 ; IWiinsch et al.|[2008l ). As far 
as infall of pristine gas clouds is concerned, our present data 
does not show any evidence for this, although the possibility 
that higher sensitivity data may show such clouds cannot be 
ruled out. 

The SBS 0335-052 system is especially unusual in that 
both galaxies in the pair are extremely metal-deficient, with 
metallicities differing by only a factor ~1.5. One possible 
explanation for their similar metallicities is that there has 
been an exchange and mixing of gas during their interac- 
tion. However, this seems unlikely, because of the faintness 
of the observed 'bridge', and also because gas exchange is 
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not believed to be very effici ent during the first encounter 
of equal-mass galaxies (|Toomre fc Toomrelll972h . It seems 
likely therefore, that the similar metallicities of these two 
galaxies are due to similar evolutionary histories. In this 
context, it is interesting to note that the SB S 0335-052 sys - 
tem is located near the edge of a large void (|Peeblesll200lh , 
and that there appears to be, in general, an indication of 
excess of metal-deficient galaxie s in voids l|Pustilnik et al.l 
l200fj ; iPustilnik fc Kniazevll2007h . 

It is also interesting to note that the most metal-poor 
galaxies known, viz., SBS 0335-052, I Zw 18 and DDO 68, 
all appear to be interacting or merger remnants. Indeed, 
in the order given above, one could take these galaxies as 
forming a 'merger sequence', (with SBS 0335-052 being 
just past the first close passage, DDO 68 being a merger 
remn ant (|Ekta et al. | l2008l; IPustilnik. Tepliakova fc Kniaze"vl 
l200Sf ) and I Zw 18 jvan Zee et al. 199Sf ) being at an inter- 
mediate phase) of extremely gas-rich galaxies. 

To summarize, our multi-resolution H I observations 
show that the SBS 0335-052 system consists of a merging 
galaxy pair. At the lowest resolution, we see a faint diffuse 
'bridge' joining the two galaxies (although the possibility 
of all of the observed emission arising from beam smear- 
ing can not be completely ruled out), while at higher res- 
olutions, we see both Hi concentrations around the stellar 
bodies as well as elongated tidal tails. The velocity fields are 
complex, and, at large-scales are dominated by the veloc- 
ity difference between the galaxies and the velocity gradient 
along the tidal tails. At still higher resolution, we see that 
in SBS 0335-052E, the ionized superbubble identified in the 
HST images by Thuan et al. (1997) is extended along one 
of the tidal tails, while there is a clump of high-density H I 
gas at the other end of the superbubble. This is consistent 
with a scenario, in which the propagating star formation 
seen in this galaxy is driven by a superbubble expanding 
into a medium with a tidally-produced density gradient. In 
contrast to SBS 0335-052E, the higher resolution maps for 
SBS 0335-052W suggest that it is undergoing a somewhat 
polar encounter. Given the similar Hi parameters of both 
galaxies, this difference in encounter geometry is possibly 
the reason for the weaker tidal disturbance and the result- 
ing large difference in the past and current star formation 
in the E and W components. 
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